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High-resolution supernova neutrino spectra represented by a simple fit 
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To study the capabilities of supernova neutrino detectors, the instantaneous spectra are often 
represented by a quasi-thermal distribution of the form fv(E) oc E a e~ f - a+ls>E ^ E '" where _B av is 
the average energy and a a numerical parameter. Based on a spherically symmetric supernova 
model with full Boltzmann neutrino transport we have, at a few representative post-bounce times, 
re-converged the models with vastly increased energy resolution to test the fit quality. For our 
examples, the spectra are well represented by such a fit in the sense that the counting rates for 
a broad range of target nuclei, sensitive to different parts of the spectrum, are reproduced very 
well. Therefore, the mean energy and root-mean-square energy of numerical spectra hold enough 
information to provide the correct a and to forecast the response of multi-channel supernova neutrino 
detection. 

PACS numbers: 14.60. Pq, 97.60.Bw 
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I. INTRODUCTION 

The neutrino signal from the next nearby core-collapse 
supernova (SN) remains the most coveted target for low- 
energy neutrino astronomy. Galactic SNe are rare, per- 
haps a few per century, and such an observation will be 
a once-in-a-lifctimc opportunity to look deeply inside a 
collapsing star and learn about its astrophysical work- 
ings as well as about neutrino properties. Several de- 
tectors worldwide are in operation that will register a 
high-statistics signal while others are in preparation or 
under discussion [1]. 

In order to assess the physics potential of various de- 
tectors or detection principles one needs the expected 
flavor-dependent SN neutrino flux spectra. In the ab- 
sence of self-consistent three-dimensional core-collapse 
simulations there is no standard SN neutrino flux model 
and moreover, the flavor-dependent flux spectra depend 
on the properties of the incompletely known neutron-star 
equation of state, on the properties of the collapsing star, 
notably its mass, and, of course, on time for any given 
case. In this situation parametric studies, taking account 
of the plausible range of predictions, are the preferred 
course of action. 

We here address one particular aspect of such studies, 
i.e. the plausible spectral shape of SN neutrino fluxes. 
On a rough level of approximation, the spectra follow a 
thermal distribution that can be described in terms of 
an effective temperature. In detail the spectra formation 
is a complicated process, different energy groups emerg- 
ing from different depths in the proto-neutron star atmo- 
sphere [2, 3], and a more refined description is necessary 
for a more detailed understanding. 

The next level of sophistication is to describe the non- 
equilibrium spectra by a three-parameter fit that al- 
lows for deviations from a strictly thermal spectrum [4]. 
One particularly simple realization are spectra where the 
energy-dependent neutrino number flux for each flavor 



has the form [3] 



f u (E) (x E a e 



a -(a+l)E/E a 



(1) 



Here, E av is the average energy and a a numerical pa- 
rameter describing the amount of spectral pinching; the 
value a = 2 corresponds to a Maxwell-Boltzmann spec- 
trum, and a = 2.30 to a Fermi-Dirac distribution with 
zero chemical potential. In general the neutrino spectra 
are fitted by 2 < a < 4 [3] with higher a > 2 indicat- 
ing stronger pinching and a < 2 meaning anti-pinching 
relative to a Maxwell-Boltzmann distribution. The three 
parameters E avi a, and the overall normalization can be 
determined, for example, if a numerical SN simulation 
provides the energy flux (luminosity) L v in some flavor, 
the average energy E av = {E v ) and some other energy 
moment, for example (i? 2 ) or sometimes (E%). 

While spectra of the form of Eq. (1) certainly pro- 
vide a reasonable overall representation, it is not obvious 
how well the spectral tails are reproduced. For study- 
ing detector responses with target nuclei with a signifi- 
cant energy threshold, it is imperative that the accuracy 
of the a-fit in the high-energy tail be checked against 
solutions of the neutrino transport equations. This is 
especially relevant for example for lead in the Halo de- 
tector [5] , argon in future large-scale liquid argon detec- 
tors [6], or subdominant detection channels on oxygen 
in water Cherenkov detectors [7] or on carbon in liquid 
scintillator detectors [8]. 

Modern numerical SN codes treat neutrino transport 
with Boltzmann solvers that are expected to produce 
physically accurate spectra. On the other hand, neu- 
trino transport is the most CPU-time consuming aspect 
of SN simulations so that in practice the energy resolu- 
tion is limited. For example, in typical simulations of the 
Garching group, 17-21 energy bins are used, correspond- 
ing to an energy resolution Ae/e ~ 0.3. Whether such 
a (seemingly) coarse zoning is adequate for modeling the 
high-energy tail of the spectrum and hence for judging 
the quality of a-fits needs to be determined at least once 
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by a resolution study. Spectra from standard resolution 
(SR) neutrino transport simulations are therefore of lit- 
tle value for assessing the quality of a-fits in the spectral 
tail, as are Monte Carlo studies [3, 4], which suffer from 
limited sampling at high energies. 

In this paper, we investigate whether a-fits and SR 
spectra from SN simulations are sufficiently accurate to 
predict detection rates that are sensitive to the high- 
energy tail. To this end, we study the spectral shape 
and its impact on the relative detector response of dif- 
ferent materials using high-resolution (HR) spectra com- 
puted for a few representative post-bounce times of a 
representative spherically symmetric simulation. These 
spectra have been obtained by re-converging models at 
a given post-bounce time after refining the energy bin- 
ning. With the HR spectra we can then test how well a 
global fit of the form of Eq. (1) reproduces the relative 
neutrino counting rates of different materials that probe 
rather different parts of the spectrum. 

We always ignore the effect of neutrino flavor con- 
versions. The ordinary Mikheyev-Smirnov-Wolfenstein 
(MSW) effect [9, 10] or self-induced flavor conver- 
sions [11] strongly modify the flavor-dependent spectra 
seen by a detector on Earth [12]. For example, the de- 
tectable v e flux will be a certain linear combination of 
the v e and u x spectra emitted at the source. These de- 
tection spectra need not be represented by a simple fit, in 
particular if self-induced conversions cause features such 
as "spectral splits." Determining the detection spectra 
in terms of the source spectra is the task of oscillation 
studies, usually performed by post-processing the fluxes 
provided by numerical models. Justifying simple repre- 
sentations of the source spectra for the purpose of flavor 
conversion studies is one prime motivation for our work. 

To this end we describe in Sec. II our numerical SN 
model and in Sec. Ill the target materials. A comparison 
between counting rates based on the numerical spectra 
and the fitted ones is performed in Sec. IV before con- 
cluding in Sec. V. 



II. NUMERICAL SUPERNOVA MODEL 

We have performed simulations of the 15M© progeni- 
tor of Woosley and Weaver [13] with the neutrino hydro- 
dynamics code Prometheus- Vertex [14], which solves 
the 0th and 1st neutrino moment equations with a vari- 
able Eddington factor closure that is provided by the 
formal solution of a simplified (model) Boltzmann equa- 
tion. The explosion of the spherically symmetric model 
was artificially initiated 500 ms after bounce. For the 
dynamical evolution of the model including the feedback 
into the hydrodynamics, we used 21 energy bins with a 
typical resolution 5e/e (ratio of zone width to zone mean 
energy) of ~ 0.29. We refer to this setup as the SR case. 

In order to accurately predict the shape of the neutrino 
spectra, particularly in the high-energy tail, we computed 
high-resolution (HR) stationary solutions of the neutrino 



transport equations for a representative time during the 
accretion phase (261 ms after bounce) and for the proto- 
neutron star cooling phase (1016 ms and 1991 ms) using 
the matter background from the dynamical simulation. 
Different from the dynamical run, the original treatment 
of [14] for the Doppler and redshift terms is used instead 
of that of [15], and the threshold values for the monochro- 
matic neutrino energy density below which a spectral ex- 
trapolation is performed has been reduced by a factor of 
100. We also transition continuously from the comoving 
frame to the lab frame once neutrino interactions cease 
in order to eliminate artifacts due to the velocity jump 
in the shock. These numerical changes accelerate con- 
vergence and help to obtain smoother and more accurate 
spectra. 

The spacing of the logarithmic energy grid was much 
finer with 42 bins and Se/e = 0.11 in the HR case. For 
better comparison, the stationary solutions for the SR 
case were re-computed with precisely the same numerics. 
Figure 1 shows the obtained HR spectra for the species 
v e , u e , and v x (with v x meaning a single kind of v ^ , v^, 
v r , and i> T ). 

In order to compare the numerical spectra (his- 
tograms) with a quasi-thermal fit of the form of Eq. (1), 
we determine the parameter a in terms of the energy mo- 
ments of the distribution. One way of fixing a is to use 
the energy moment of order k and notice that 



(E k ;) 
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i 



(Eu). 



(2) 



We have checked that for our cases the implied value of 
a does not strongly depend on k, which already implies 
that this type of fit is a reasonable representation. In the 
following we will use k = 2, i.e., the relation 
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We show examples of the fits as continuous colored curves 
alongside the HR spectra (step functions) on a logarith- 
mic and on a linear scale in Figs. 1 and 2, respectively. 
For comparison we also show the SR spectra and their 
correspondent fits as thin dashed black lines in Fig. 1, and 
find that visually the fits work very well in both cases. 
Values for the fit parameter a, the luminosity L„, and the 
neutrino mean energy {E„) for the three snapshots and 
for the different neutrino species v e , v e , and v x are listed 
in Table I both for the HR and the SR case. We note 
that despite the lower energy resolution, the SR spectra 
from the dynamical simulation typically give very simi- 
lar luminosities and mean energies, which are consistent 
with the HR case within ~ 5% or less. The fit parameter 
a agrees similarly well except for v e in the post-explosion 
phase with a deviation of up to 9%. However, this is due 
to a very steep dependence of a on the second energy 
moment, (E„), when the energy ratio in Eq. (3) is close 
to unity, and does not affect the fit curves appreciably. 
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FIG. 1: Spectra for species v e (left column), (middle column) and v x (right column) for post-bounce times of 261 (top row), 
1016 (middle row), and 1991 ms (bottom row). The data from the Prometheus- Vertex spectra are shown as step functions, 
and the continuous curves are quasi-thermal fits according to Eq. (1). Thick magenta, turquoise, and green lines are used for 
the different neutrino species (y e , i> e and v x , respectively) in the HR case, while the SR spectra and fits are shown as thin 
dashed black lines. 



TABLE I: Fit parameters a, luminosities L v , and neutrino mean energies (E v ) for different neutrino species and post-bounce 
times for HR and SR spectra. 



t 


species 


Qhr 


c*sr 




ii/,SR 


(-Bi/)hr 


(^)SR 


[ms] 








[10 52 ergs" 1 ] 


[10 52 ergs" 1 ] 


[MeV] 


[MeV] 


261 


v e 


2.65 


2.73 


2.80 


2.68 


13.05 


12.87 


261 


v e 


3.13 


3.23 


2.79 


2.71 


15.23 


15.05 


261 


v x 


2.42 


2.48 


1.58 


1.63 


14.48 


14.48 


1016 


Ve 


2.90 


3.14 


0.62 


0.62 


10.14 


10.28 


1016 


V e 


2.78 


2.88 


0.67 


0.70 


12.69 


12.87 


1016 


Vx 


2.39 


2.52 


0.75 


0.77 


12.89 


13.00 


1991 


V e 


2.92 


3.17 


0.48 


0.47 


10.01 


10.01 


1991 


V e 


2.61 


2.67 


0.44 


0.44 


12.28 


12.36 


1991 


v x 


2.34 


2.49 


0.54 


0.55 


12.31 


12.39 



III. DETECTION CROSS SECTIONS (IBD) reaction v e +p -> n + e+. The cross section varies 

roughly as E% so that the overall event rate is well re- 
in existing detectors, SN neutrinos are primarily mea- produced, if the SN neutrino spectrum is approximated 
sured in the D e channel by virtue of the inverse beta decay 



4 



CD 




10 20 30 40 50 

E (MeV) 



FIG. 2: Spectra for species v^, i/ e and v x . Histograms are 
the Prometheus- Vertex HR spectra of our model at time 
1016 ms. The curves are quasi-thermal fits according to 
Eq. (1): dashed for v & , continuous for £> e , and dot-dashed 
for v x . 
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FIG. 3: Neutrino cross sections on some representative target 
nuclei (see text for details). 



by a quasi-thermal spectrum where the second energy 
moment is used to determine a. 

However, other detection channels are also of interest. 
One is elastic scattering on electrons where the cross sec- 
tion varies roughly linearly with E v . Of greater interest 
for our purpose of testing the goodness of the a-fit are 
reactions with a higher energy threshold or steeper en- 
ergy variation than IBD. One example is the recently 
built HALO detector in SNOLAB where 79 tons of lead 
are used as target material. Other examples are reac- 
tions on argon in future liquid argon detectors that are 
primarily considered in the context of long-baseline neu- 
trino oscillation studies and reactions on oxygen in water 
Cherenkov detectors and on carbon in liquid scintillator 
detectors. 

In our study we consider several selected reactions that 



span a broad range of spectral responses to neutrino 
fluxes and show the energy-dependent cross sections in 
Fig. 3. The reaction sensitive to the lowest range of en- 
ergies is elastic scattering on electrons via both charged- 
current (CC) and neutral-current (NC) for v e and v e and 
via NC interaction for v x [16]. 

In nuclei, CC interactions proceed via v e + (N, Z) — > 
(N-l,Z+l) + e~ and D e + (N,Z) -> {N + l, Z-l) + e+, 
respectively. The v e interaction is typically suppressed 
at a given energy with respect to the v e interaction due 
to Pauli blocking. Neutral current interactions may also 
produce observable signals via ejected nucleons or de- 
excitation photons. As one example we consider NC 
interactions with 16 [7, 17-19], which is relevant for 
water-based detectors. 

Interactions with heavier nuclei, such as lead, may 
yield quite high rates of both CC and NC interactions. 
Observable signatures include leptons and ejected nucle- 
ons. Single and multiple neutron ejections are possi- 
ble. The relevant interactions for lead-based detectors 
are v e + A Pb -> e" + ^Bi* and v x + A Pb -> v x + A Pb* . 
For both CC and NC cases, the resulting nuclei de-excite 
via neutron emission. Antincutrino CC interactions are 
strongly suppressed. Although natural lead contains iso- 
topes other than 208 Pb, the neutron cross section for 
A = 208 should be similar to other components. For 
the Pb cross sections see Refs. [20-23]. 

Liquid argon detectors will have excellent sensitivity to 
v e via the CC interaction j/ e + 40 Ar -> e" + 40 K*. This is 
an interaction for which the de-excitation 7s from 40 K* 
can be observed. The reaction v e + 40 Ar — > e + + 40 C1* 
will also occur and can be tagged via the pattern of 7s. 
NC excitations are possible, although little information 
is available in the literature because of the low-energy 
nuclear recoils that make it practically difficult to detect. 
Energy thresholds as low as few MeV may be possible. 
Although the chances for detecting such a reaction are 
rather low, we consider NC coherent scattering anyway as 
an example of a reaction scanning the low-energy region 
of the neutrino spectra. The cross sections are reported 
in Refs. [6, 24-27]. 

For the different channels that wc consider, we use 
the cross sections as in [1]. They have been used in the 
SNOWGLOBES software package [28]. 



IV. COUNTING RATES 

The detector counting rate as a function of incoming 
neutrino energy is proportional to f v (E)cr(E). Based on 
the numerical neutrino spectra for the 1016 ms model, 
we show this quantity in Fig. 4 for our suite of interac- 
tion processes, normalized to the maximum rate of each 
process. This plot demonstrates that the different detec- 
tion processes probe vastly different parts of the neutrino 
spectra. 

We now compare the total counting rates for each of 
our reactions using as a neutrino spectrum either the 
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FIG. 4: Detection rate as a function of neutrino energy for 
different reactions, based on our 1016 ms model. Each curve 
is normalized to its maximum. Top: Anti-neutrinos. Bottom: 
Neutrinos. 



numerical output directly (rate N num ) compared with 
the rate inferred if the spectrum is represented by our 
fit function (rate A^ t )- In Table II, we report the ra- 
tio ./Vfit /-/Vnum for neutrinos (top part) and antineutrinos 
(bottom part). The agreement between N mm and Nat is 
very good and typically the error is only a few percent. 
Even for the oxygen reaction, which probes the highest- 
energy tail of all examples, the difference is only some 
10-20%. Assuming complete flavor conversions due to 
neutrino oscillations (v e — > v x , and the same for antineu- 
trinos), we tested that the ratio Nfi t /N num differs from 
the ones reported in Table II only on the level of tenths of 
a percent (results not shown here) . Such a result further 
confirms the good quality of the alpha-fit. 



V. CONCLUSIONS 

In this paper, we have investigated whether the sim- 
ple analytic fit of Eq. (1), based on the first two energy 



moments of the numerical distribution, is sufficiently ac- 

TABLE II: Ratio of counting rate based on our spectral fit 
and the direct numerical HR spectrum, Nat/N nunl . 



t[ms] 


z/ e 40 Ar 


^ e 208 pb 


f e e 


v x e 






CC 


CC 






NC 


261 


1.01 


1.02 


1.00 


1.00 


1.01 


1016 


1.02 


1.06 


1.01 


1.00 


0.91 


1991 


1.02 


1.07 


1.01 


1.00 


0.90 




t\ms] 


P e 40 Ar 


^ e 208 Pb 


v e e 


v x e 


* e 16 




NC 


NC 






NC 


261 


1.00 


1.02 


1.00 


1.00 


1.05 


1016 


1.01 


1.03 


1.00 


1.00 


1.16 


1991 


1.00 


1.03 


1.00 


1.00 


1.18 



curate to predict the detector response for target nu- 
clei with relatively high threshold energies. In order to 
answer this question reliably, we have produced high- 
resolution (HR) neutrino spectra at several postbouncc 
times of a spherically symmetric SN simulation. In ad- 
dition, these spectra have allowed us to check for the 
first time whether the coarser standard resolution (SR) 
of typical supernova models (which is dictated by strong 
CPU-time limitations) already allows robust signal pre- 
dictions. 

We have compared SR and HR spectra and fitted both 
with the analytic fit of Eq. (1), based on the first two en- 
ergy moments of the numerical distribution. These fitted 
spectra account well for the detection rates in SN neu- 
trino detectors with vastly different target nuclei. Wc 
conclude that for the purpose of signal forecast in differ- 
ent detectors reasonably good accuracy can be achieved 
by considering the lowest two energy moments of numer- 
ically produced spectra. Moreover, we have verified that 
a modest standard energy resolution as affordable in typ- 
ical supernova simulations already gives energy moments 
surprisingly close to the HR case. This insight simplifies 
parametric studies of detection forecasts for the neutrino 
signal from the next nearby SN because it is sufficient to 
represent the non-equilibrium spectra by a simple three- 
parameter quasi-thermal fit. 
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